	/*Reproduce main model in Stata*/
		use nonimp,replace

	set more off

	/*Drop Luxemburg*/
	drop if zreplacementrate ==.
	/*Drop all cases with missing micro data*/
	drop if rexvote == .
	drop if male == .
	drop if age1 == .
	drop if age2 == .
	drop if age4 == .

	drop if myeduc == .
	drop if myoccup == .
	drop if zlrs == .
	drop if zsatisdmo == .
	drop if euschlecht == .

	/*generate dummies for countries*/
	desmat sortcountry ,full

	/*generate interaction effects*/
	
	gen intasylarbeit=zasylumseekers * zsur
	gen intzreplacezsur = zsur*zreplacementrate
	gen intzreplacezasyl = zasylumseekers*zreplacementrate

	/*center manifesto variables*/
	replace salienzmean = salienzmean - 3.84568
	replace rvar = rvar - 21.75423
	gen insalienzrvar = salienzmean * rvar

	/*Estimate model - go for a large cup of coffee!*/
	xtmelogit rexvote male age* mye1 mye2 farmerown -unemployed zlrs euschlecht zsatisdmo disp lfed1 zasylumseekers zsur intasylarbeit zreplacementrate intzreplacezsur intzreplacezasyl rmax salienzmean rvar insalienzrvar _x*,nocons || kontext: , cov(identity)

	est save fullmodel, relace
